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Abstract - We introduce a new method for detection of long-range cross-correlations and mul- 
tifractality - multifractal height cross-correlation analysis (MF-HXA) - based on scaling of qth 
order covariances. MF-HXA is a bivariate generalization of the height-height correlation analysis 
of Barabasi & Vicsek [Barabasi, A.L., Vicsek, T.: Multifractality of self- a ffine fractals, Physical 
Review A 44(4), 1991]. The method can be used to analyze long-range cross-correlations and 
multifractality between two simultaneously recorded series. We illustrate a power of the method 
on both simulated and real- world time series. 
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i_i The research of long-range dependence and multifrac- 
tality has been growing significantly in recent years with 
^^application to a wide range of disciplines [I 10 . Recently, 
^y^the examination of long-range cross-correlations has be- 
f^come of interest as it provides additional information 
j" about the examined processes. Carbone [II] generalized 
CO the detrending moving average (DMA) method for higher 
dimensions. Podobnik & Stanley 12 adjusted the de- 
trended fluctuation analysis for two time series and intro- 
duced the detrended cross-correlation analysis (DCCA). 
. . Zhou [13] further generalized the method and introduced 
!>the multifractal detrended cross-correlation analysis (MF- 
k>DXA). Jiang & Zhou [14] then implemented moving aver- 
^ age filtering to MF-DXA algorithm creating MF-X-DMA. 
Cdln this paper, we introduce two new methods for an 
analysis of long-range cross-correlations - the multifractal 
height cross-correlation analysis (MF-HXA) and its spe- 
cial case of the height cross-correlation analysis (HXA). 

To analyze long-range cross-correlations, we generalize 
the q-th order height-height correlation function for two 
simultaneously recorded series. Let us consider two series 
X t and Y t with time resolution v and t — ^,2^, ...,z/|_— J, 
where [J is a lower integer sign. For better legibility, 
we denote T* = ^[77 J? which varies with v, and we 
write the r-order difference as A T X t = X t+T — X t and 
A T X t Yt = A T X t A T Y t . Height-height covariance function 
is then defined as 



T*/v 



K xy , q (r) = —J2 \ A rX t Y t \i = (\A T X t Y t \i) (1) 



t=i 



where time interval r generally ranges between v = 
T~min, • • • j T max . Scaling relationship between K xy ^ q (r) and 
the generalized bivariate Hurst exponent H xy (q) is ob- 
tained as 



.^Wocr^W. (2) 

For q = 2, the method can be used for the detection 
of long-range cross-correlations solely and we call it the 
height cross-correlation analysis (HXA). Obviously, MF- 
HXA reduces to the height-height correlation analysis of 
Barabasi et al. [l5] for X t = Y t . Note that it makes sense 
to analyze the scaling according to Eq. [2] only for de- 
trended series X t and Y t and only for q > pi. A type 
of detrending can generally take various forms - polyno- 
mial, moving averages and other filtering methods - and 
is applied for each time resolution v separately. 

The bivariate Hurst exponent < H xy (2) < 1 has simi- 
lar properties and interpretation as a univariate Hurst ex- 
ponent. For H xy (2) > 0.5, the series are cross-persistent so 
that a positive (a negative) value of AX t AY t is more sta- 
tistically probable to be followed by another positive (neg- 
ative) value of AXt+iAlt+i- Conversely for H xy (2) < 0.5, 
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the series are cross-antipersistent so that a positive (a neg- 
ative) value of AX t AY t is more statistically probable to be 
followed by a negative (a positive) value of AX t +iAY t +i. 
Note that even two pairwise uncorrelated processes can be 
cross-persistentQ 

The expected values of the bivariate Hurst exponents 
have been partly discussed in 
that 



12 -14 . It has been shown 



H xy(q) 



H x (q)+H y (q) 



(3) 



for all q > for pairwise uncorrelated and correlated 
processes. We present some new insights into this rela- 
tion. To better understand the behavior of the bivariate 
Hurst exponent, we use a standard multifractal formalism 
[16] . Consider processes X t and Y t are multifractal with 
generalized Hurst exponents H x (q) and H y (q) so that 



(\A T X t \ q ) (xr qHM 
(\A T Y t \ q ) (xr qH ^ q \ 



(4) 
(5) 



In the same way, we can write the joint scaling of two 
series (compare with Eq. [I]) as 



\A T X t Y t \i) oc T qH ^ q \ 



(6) 



Using the definition of covariance, the left part of Eq. [6] 
can be rewritten as 

{\A T X t Y t \i) = {\A T X t \i){\A T X t \i) + cov(\A r X t \i ,|A T Y t |§). 

(7) 

From Eqs. [4]and[5j the first part of the right-hand side 
of Eq. [7] implies 

q q H x (q) + H v (q) 

(|A T X t |«)(|A T F t |f)ocr« ^~ (8) 

which corresponds to Eq. [3j Therefore, the crucial part 
of long-range cross-correlations and multifractality is the 
scaling of covariances between \A r X t \i and |A T >t|2 with 
varying r. Consider now a scaling exponent a(q) and a 
scaling relationship 



C0V (|A T X t |3,|A T F t |3)cXT 



qa(q) 



(9) 



This leads us to three simple implications. If covariances 
do not scale with r, then Eq. [3] holds. If the covariances 
scale with r, the other two are as follows: 



1 For example, let us have pairwise uncorrelated processes Xt 
and Yt following fractional Gaussian noise with H x (2) = 0.9 and 
H y (2) = 0.7 and thus H xy {2) = 0.8. Even though the two pro- 
cesses are independently generated (and thus uncorrelated), they 
are cross-persistent. If Xt > and Yt > 0, then it is statistically 
more probable (based on persistence of the separate processes) that 
also Xt+i > and lt+i > than otherwise. Therefore, if the pro- 
cesses moved together in period f, it is statistically more likely that 
they will move together in period t + 1 as well (and vice versa), i.e. 
the processes are cross-persistent. 



a(q) 



H x (q) + H y (q) 



H xy (q) 



H x (q) + H y (q) 



(10) 

We show that these relationships are indeed true for ar- 
tificially generated processes later in the text. Therefore, 
we need to distinguish between two types of long-range 
cross-correlations: (i) long-range cross-correlations caused 
by long-range dependence of the separate series, and (ii) 
long-range cross-correlations caused by scaling of covari- 
ances between |A r X t |2 and \A T Y t \z . 

In order to test validity of the method, we present re- 
sults for several artificial series. In the analysis, we ap- 
ply MF-HXA with changing r max = 5, . . . , 100 and fixed 
Tmin = 1. In turn, we obtain the 99% jackknife confidence 
intervals under an assumption of a normally distributed 
Hurst exponent with an unknown variance. The estimated 
Hurst exponent is then taken as a mean of the exponents 
based on the various r max . This way, we can comment on 
the results with statistical power 5 . In the procedure, we 
apply filtering of a constant trend. We now turn to the 
artificial processes. 

First, we start with the Mandelbrot's binomial multi- 
fractal (MBM) measures 



17 



18 . Let m > 0, mi > 



and mo + mi = 1 and let us work on interval [0,1]. In 
the first stage, the mass of 1 is divided into two subin- 
tervals [0,|] and [|,1], when there is the mass mo in the 
first subinterval and the mass mi in the second one. In 
the following stage, each subinterval is again halved and 
its mass is divided between the smaller subintervals in a 
ratio mo : mi. After k stages, we obtain a series of 2 k 
values. Note that the values are deterministic as there 
is no noise added in the simplest version of the method. 
For an interval [z,z + 2 _fe ], the value \i has a value of 
n[z, z + 2~ k ] = mj^m^ 1 , where cpo and cpi stand for the 
relative frequencies of numbers and 1 in a binary devel- 
opment of 2 k z, respectively. We construct two series with 
mo = 0.3, 0.4 and k = 16. Results are presented in Fig. 
|TJl, showing that the bivariate Hurst exponent H xy (q) does 
not deviate significantly from the average value of H x (q) 
and H y (q) even though the analyzed series are strongly 
correlated. 

Second, we apply MF-HXA on ARFIMA processes with 
correlated noise terms. ARFIMA(0, 0) process is defined 
as x t = YliLi Q>i(d)xt-i + St where < d < 0.5 is a free 
parameter, related to Hurst exponent as H = d + 0.5, and 
a,i(d) = dT(i - d)/(r(l - d)r(l + z)). We simulate long- 
range dependent series of length 10 4 . To describe influence 
of the correlations on MF-HXA estimates, we generate se- 
ries with correlated noise e t ~ 7V(0, 1) and five cases are 
investigated - correlation coefficients for the noise terms 
are set to 1, 0.5, 0, -0.5 and -1. The results are shown 
in Fig. [TJ> f. The estimates of H xy (q) are not signifi- 
cantly different from Hx ^~^ H y^ for any q or any corre- 
lation coefficient value. This result is in hand with the 
results shown in [l4] - pairwise correlations have no effect 
on the H xy estimation. 

Third, we analyze the behavior of two-component 
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ARFIMA processes [19]. For parameters d\ and c?2, the 
two-component ARFIMA(di,c?2) processes X t and Y t are 
described by the following set of equations: 

X t = [Wx t + {l-W)yt]+e t 
Y t = [(1 - W)x t + Wy t ] + vt 

x t = Y^Li a ii. d i) x t-i 
yt = ^2iLi a i( d 2)Y t -i 

Here, W is a free parameter (0.5 < W < 1) controlling 
a strength of coupling between X t and Y t , and s t ,vt ~ 
iV(0, 1) are noise terms. Note that for W = 1, we obtain 
two decoupled ARFIMA processes, whereas for W < 1, 
the two processes have long memory of the process itself 
as well as of the other one. In our simulations, we con- 
sider di = c?2 = 0.3 with W = 0.5,0.75 (practically, the 
case W = 1 has been investigated in the previous para- 
graph). The results are shown in Fig. [T^h. For both 
W = 0.75 and W = 0.5, we notice deviations of H xy (q) 
from ILM+EmM starting already at q = 0.1. The devi- 
ations are statistically insignificant for lower moments q 
(due to rather short series, T = 10 4 ), but become statis- 
tically significant for higher moments (for q > 1.3 when 
W = 0.5 and for q > 2.5 when W = 0.75). The effect gets 
stronger with lower W. Indeed, these are expected results 
as the construction of the two-component ARFIMA mixes 
the long memory of the separate processes together. 

In Fig. [2| we present the results based on separation 
in Eq. [7| i.e. scaling of separate processes and scaling 
of covariances of \A r X t \z and |A r l^|2. For illustrational 
purposes, we show only the case q = 2. For MBM (Fig. 
[2|l), the scaling of covariances is slightly lower than the av- 
erage of Hurst exponents, yet remains well between them. 
This is reflected in the fact that the estimated H xy (2) is 
not equal to the average of estimated H x (2) and H y (2) 
but is rather close to the lower confidence interval (Fig. 
|TJl) , yet the deviation is still insignificant. In Fig. |2|d, four 
cases of correlated ARFIMA processes are illustrated. All 
four processes show a(2) w 0.7, which perfectly fits the 
expectations. We can see that the covariances are higher 
for highly correlated series than the less correlated series, 
but the scaling relation remains the same for all. The case 
of uncorrelated ARFIMA processes exhibits no scaling of 
covariances (as these vary around zero) and is thus not 
shown. In Figs. [2j;,d, the two-component ARFIMA pro- 
cesses are illustrated. Here, the difference between scaling 
of covariances and the pair of K x ^(t) and K y ^(r) is re- 
markable for both W — 0.75 and W = 0.5. The scaling of 
covariances is expectedly stronger for W = 0.5. These re- 
sults perfectly support the calculations presented in Eqs. 

ED uni 

To show potential use of the method, we study differ- 
ent real-world financial series, which we consider the out- 
puts of the complex systems - daily volatility and vol- 
ume series of NASDAQ and S&P500 stock indices (fi- 
nance. yahoo. com database), and daily returns and volatil- 
ity of spot and futures prices of WTI Crude Oil (NYMEX 



Commodities database). Even though the real- world se- 
ries are of the same length order as the simulated pro- 
cesses, which scale even up to r = 100 and q — 10, K T (q) 
usually does not scale for r > 20 and q > 3 for daily finan- 
cial data [5]. Also, we apply linear filtering according to 
(4]. The generalized Hurst exponents are then estimated 
by varying T max between 5 and 20 for 0.1 < q < 3. 

For the stock indices, we analyze the series of vol- 
ume and volatility for the longest available datasets - 
from 11.10.1984 to 26.4.2011 for NASDAQ (6,693 ob- 
servations) and from 3.1.1950 to 26.4.2011 for S&P500 
(15,428 observations). We take absolute returns, defined 
as I log P t — log Pt-i \ where P t is a stock index closing price, 
as a measure of volatility. Volume series are transformed 
as a relative deviation from a moving average of traded 
volume in approximately past two trading years (500 ob- 
servations) to control for an exponential increase of the 
traded volume in past decades (Fig. [3^i,c). The estimated 
generalized Hurst exponents are shown in Fig. [3]3,d. For 
both stock indices, the trading volume and volatility are 
strongly persistent as well as cross-persistent. Neverthe- 
less, the bivariate Hurst exponent H xy (q) does not differ 
significantly from the average of H x (q) and H y (q), i.e. the 
cross-persistence of the series is mainly due to the persis- 
tence of the separate processes and the fact that the pro- 
cesses are correlated (Fig. [4^,f). The scaling of K xyjq (r) 
is very stable up to r = 20 and for all examined qs (Fig. 
[4^i, b). The results are in hand with |12 who found weaker 
persistence of the process of traded volume. However, the 
definitions of traded volume differ from our study. 

For the WTI spot and futures prices, we cover a period 
from 2.1.1986 to 26.4.2011 (6,348 observations) and ana- 
lyze the logarithmic returns r t = \og(P t ) — log(Pt-i) and 
the volatility again in the form of absolute returns. The re- 
sults are shown in Fig. |3^,f. For both returns and volatil- 
ity, the estimates of the generalized Hurst exponents prac- 
tically overlap for all q. On one hand, the returns show no 
signs of long-range correlations or cross-correlations. On 
the other hand, the volatility of separate processes show 
strong persistence as well as cross-persistence. Moreover, 
the generalized Hurst exponents vary only slightly with 
q and are not even monotonically declining as expected 
for multifractal processes, suggesting that the processes of 
volatility are monofractal. Yet again, the cross-persistence 
of the series is mainly due to the persistence of the sepa- 
rate processes and high correlation between the processes 
(Fig. 4g,h) as H xy (q) does not significantly deviate from 



9 . The scaling of K xyA {r) shows different be- 
havior for returns and volatilities. As for volatility, the 
scaling is very stable up to r = 20 and q = 3. On contrary, 
the scaling for returns becomes less stable with growing q 
(Fig. §,d). 

In conclusion, we introduce the new method for an anal- 
ysis of long-range cross-correlations and multifractality 
- the multifractal height cross-correlation analysis. The 
scaling of covariances of the absolute values of the series 
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gives additional information about dynamics of two simul- [19] Podobnik B., Horvatic D., Lam Ng A., Stanley H. 
taneously recorded series and can cause divergence of the and Ivanov P., Physica A , 387 (2008) 3954. 

bivariate Hurst exponent from the average of the sepa- 
rate univariate Hurst exponents. A utility of the method 
has been shown on several artificial series as well as the 
real- world time series. We argue that even though major- 
ity of the analyzed series are cross-persistent, such cross- 
persistence is mainly caused by persistence of the sepa- 
rate processes and the fact that the series are correlated. 
The scaling of covariances of the absolute values of the 
examined processes is with good agreement with this re- 
sult. A larger study comparing bias and efficiency of MF- 
HXA compared to the other methods analyzing long-range 
cross-correlations (MF-X-DFA and MF-X-DMA) shall fol- 
low. 
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Fig. 1: (a) Binomial multifractal measures. Generalized Hurst exponents (y-axis) dependent on moments q £ [0.1,10] (x-axis) with step 
of 0.1. H x (q) for MBM with mo = 0.3 (bold black line) varies only weakly with q compared to H y (q) for MBM with mo = 0.4 (bold 
dashed black line). H xy (q) (bold red line) is not significantly different from the average of H x and H y (dotted black line) for all q (the 
99% jackknife confidence intervals around H xy (q) in gray), (b) - (f) ARFIMA(0,d,0) processes with correlated noise. Generalized Hurst 
exponents (y-axis) dependent on moments q £ [0.1,10] (x-axis) with step of 0.1. H x (q) for ARFIMA(0,d,0) with H = 0.8 (bold black 
line) and H y (q) for ARFIMA(0, <i, 0) with H = 0.6 (bold dashed black line). The rest of the notation and parameters setting holds from 



0.5, 



0, p 



-0.5 



(a). Figs, (b) - (f) show ARFIMA(0, d, 0) processes with correlated noise with correlation coefficients p = 1, p 
and p = —1, respectively. There is no significant deviation of H xy (q) from Hx ^~^ H y^ f or a n q anc j f or a \\ examined correlations, (g) - 
(h) Two component ARFIMA(d\, cfoj processes. Generalized Hurst exponents (y-axis) dependent on moments q £ [0.1,10] (x-axis) with 
step of 0.1. Here, we use two component ARFIMA(di , cfo) processes with d\ = cfo = 0.3 and varying W. For W = 0.75 (g), H xy {q) is 
significantly higher than Hx ^~^ H y^ f or higher moments (q > 2.5). For W = 0.5 (h), the deviation of H xy (q) from Hx ^~^ H y^ j s higher 
than for case (g) and the statistically significant deviation from the average starts at lower moments (q > 1.3). 
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Fig. 2: Scaling of K x ^(j), K y ^(j) and covariances between \A T Xt\ and |A T Yt|. The scaling is shown for MBM (a), correlated ARFIMA 
processes (b) and two-component ARFIMA processes (c,d). For MBM, we observe slight divergence of a (2) from the average of H x (2) 
and H y (2), which remains insignificant (see Fig. ^i). Correlated ARFIMA processes show practically perfect fit the expected a(2) of 0.7 
(ARFIMA processes with d = 0.1 and d = 0.3). Two-component ARFIMA processes exhibit remarkable deviation of a from the average 
of Hurst exponents. Note that the fits (dashed black lines) and slopes are estimated on the whole sample from r m i n = 1 to r m ax = 20. 
The results are in hand with expectations based on Eqs. [4]-|10|and in agreement with Fig. ^ 
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Fig. 3: (a), (b) NASDAQ, (a) Evolution of logarithmic prices (red, left y-axis) and transformed traded volumes (black, right y-axis) in 
between 11.10.1984 and 26.4.2011. (b) Generalized Hurst exponents (y-axis) dependent on moments q £ [0.1,3] (x-axis) with step of 0.1. 
H x (q) for NASDAQ volatility (bold black line) and H y (q) for NASDAQ traded volume (bold dashed black line) both vary with q while 
stronger variation is present for volume. H xy (q) (bold red line) is not statistically different from the average (dotted line) of H x (q) and 
H y (q) for any q. (c), (d) S&P 500. The time period covered ranges from 3.1.1950 to 26.4.2011. Same notation and estimation parameters 
setting hold here. H xy again does not differ from the average of the univariate Hurst exponents, (e) WTI crude oil spot and futures 
prices returns. Same notation holds, H x represents the dynamics of spot returns and H y for futures returns. There is again no significant 
deviation of H xy (q) from Hx ^~^ H v^ m (f) WTI crude oil spot and futures prices volatility. The notation holds. Generalized Hurst 
exponents practically overlay for both series as well as for the joint dynamics. 




Fig. 4: (a) - (d) Scaling of K xy , q for NASDAQ, SPX and WTI. The scaling functions are presented by black lines, where the cases of 
q = 1,2,3 are in bold. For these three cases, the best fits are illustrated (dashed lines). The scaling is very stable for NASDAQ, SPX 
and WTI volatility up to r = 20 for all examined q so that the fits are almost undistinguishable from the scaling functions. For WTI 
returns, the scaling is less stable with increasing q. For the analyzed series, it implies that scaling is better for higher values of Hurst 
exponents, (e) - (h) Scaling of K x ^ q , K y ^ q and covariances between \A T Xt\ and \A T Yt\ for NASDAQ, S&P500 and WTI. Best linear 
fits are represented by dashed lines and estimated slopes are noted. For illustrational purposes, we show only the case q = 2. The scaling 
exponents a (2) are approximately equal to the average of estimated Hurst exponents. This implies that eventual cross-persistence (for 
cases of NASDAQ, S&P500 and WTI volatility) is majorly caused by persistence of the separate processes and the fact that the processes 
are pairwise correlated. Note that the differences of estimates from Fig [3] are caused by the fact that here, we estimate the exponents for 
r between r m i n = 1 and r m ax = 20, while for Fig. [3] we use the jackknife estimates. 
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